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Abstract 



We propose a theory for the observed tricriticality in the orientational phase 
diagram of Si(113) misoriented towards [001]. The systems seems to be at or 
close to a very special point for long range interactions. 
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A basic question to all surface studies is the stability of a surface against various processes, 
e.g., thermal fluctuations, defects, step formation etc. This question assumes importance 
especially if the surface is cut, not in a perfect crystallographic direction, but with a slight 
miscut angle (called misorientation). So far, miscut surfaces have been found to be stable 
with steps or to reorganize to more complex surface structures [|TJ. An exception to this 
general rule is the recent observation of a tricritical point for a Si(113) surface with a small 
miscut angle towards [001], point at a temperature T t = Y12ZK. A phase coexistence 
of a stepped surface with a (113) facet, for T < T t was observed in Ref 0|§. For T > 
T t , the surface goes continuously to the facet as the misorientation is decreased. In plain 
English, the crystal can be cut at any small angle as one wishes so long T > T t but, for 
any T < T t , there is a minimum miscut angle below which the crystal surface cannot be cut 
( in thermal equilibrium). The phase coexistence of a reorganized facet and steps can be 
understood on the basis of two separate free energy curves |lj but, a tricritial point demands 
a more subtle treatment. This discovery of a tricritical point in a two dimensional system is 
extremely important because it can serve as a fertile ground for recently developed statistical 
mechanical theories [§|-|9|]. Our aim is to describe the universal aspects of this tricritical point. 

The steps run in one direction ("z" axis) without backtracking (no overhang), so that the 
surface can be characterized by the density of steps, p (number of steps/transverse length) 
Q. The misorientation is related to p. See Fig 1. In the experiment of Ref. @|§, the 
coexistence curve [p ~ (T t — Tf ], was found to have a zero slope at the tricritical point, 
with (3 = 1/2, showing utter disrespect to the extant mean field theory ||[| that predicts 

The phenomenology of the transition can be discussed in terms of the Legendre transform 
of the free energy per unit transverse length, F(p, T) 

F(ji,T;p) = -w + F(p,T), (1) 

where p is the chemical potential for steps. The thermodynamic value of p comes from 
dJ- '/dp = 0, or equivalently dF/dp = p. Taking f as the free energy of an isolated step, 
Eq. (1) can be expressed as 

T; p) = (f - p)p + F int {p, T), (2) 

where Fi nt (p,T) is the free energy contribution from interactions among the steps. For 
noninteracting steps, Fi nt = 0. Therefore, a first order transition at p = fo takes a facet 
(p = 0) to a fully stepped surface for p > f . For purely repulsive ( "fermionic" ) steps, F int ~ 
p 3 leading to the continuous Pokrovsky-Talapov (or "3/2" order) [|l^|ll],||,[l!| transition with 
P H fo^f 1 l 1 ^ 2 - This has been the rule for almost all systems until the tricriticality in Si(113) 
was discovered. 

The occurrence of a phase separation, as noted in ref ||, suggests the existence of at- 
tractive interactions among the steps. The steps are generally taken to be nonintersecting. 
In addition, they are expected to have dipolar or elastic long range r~ 2 interaction. A 
long range (r~ 2 ) attraction for the non-intersecting steps, in a mean field ("Hartree Fock") 
analysis, gives F int ~ (T — T t )p 3 + ap 4 , predicting a tricritical point with a linear phase 
boundary p[. With an attractive r~ 2 interaction, this mean field form of Fi nt agrees, in the 



domain of overlap of parameters, with the exact Bethe ansatz solution [13]. This rules out 
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the possibility of fluctuations causing a zero slope phase boundary in this model with long 
range attration. 

We therefore consider a completely different scenario. Our proposal is that the tricrit- 
icality occurs when the steps have a short range attraction. The physical picture we have 
in mind has strong resemblance to a phase separating polymer solution. In fact taking the 
steps as directed polymers (DP), the stepped face can be thought of as a DP solution. Our 
proposal is that in the high temperature phase, the steps are repulsive and the phase is 
dominated by the entropic interaction. As the tricritical point is approached, the attractive 
part starts playing a role. The steps start colliding and the average separation between 
collisions determine the correlation length. When this length becomes comparable to the 
separation of the steps the tricritical domain is reached. The phase separation takes place 
in this regime before any bound state can form. 

To study the phase separation and the coexistence curve, we use a canonical ensemble 
approach, and use the analogy with polymer solution fT7|] . The phase boundary is identified 
by equating the "osmotic" pressure of the two coexisting phases. The "osmotic" pressure in 
this context, would mean the excess pressure generated by the addition of one more step and 
is obtained from II = p 2 d/dp(F/p), where F(p) is the free energy as a function of density 
p. Since one phase is a flat surface with zero density, its osmotic pressure is zero. The 
coexistence curve is therefore obtained from 

d (p^F) /dp = 0. (3) 

It is the interaction that determines the shape of the coexistence curve. We, therefore, 
consider two different possibilities: tricriticality with (i) short range interactions, and (ii) 
long range interactions. It seems that the latter holds the key. 

The general approach is to start from a mean field or effective free energy. Renormaliza- 
tion group (RG) approach is then used to incorporate the effects of fluctuations. The RG 
f3 functions tell us the effective couplings as the length scales are changed. Integrating the 
RG equations, one can then obtain the renormalized interactions or coupling constants for 
the relevant length scale ~ p _1 . These renormalized coupling constants can then be used in 
the mean field free energy to get F for Eq. [5| This is justified because we are interested not 
in the details of the tricritical behavior but rather in the phase boundary where all length 
scales remain finite. 

To write the Hamiltonian, we note that short range attractions in DP's lead to bound 
states for two isolated steps. So far as the binding transition is concerned in low dimen- 
sions, the universal critical behaviour is independent of any further details like noncrossing 
condition of the steps. This is known from exact renormalization group analysis and simple 



quantum mechanical calculations |14l[T6],|T5| . Treating the steps as structureless wandering 



lines, the Hamiltonian in a continuum formulation is taken as 



H — dz 




+ V 2 J2 6 ( r ij( Z )) 

i>j 



+ H int , (4a) 



where Ti(z) is the d dimensional transverse position of the ith step at a coordinate z measured 
along the step from one end, (2) = Vj{z) — Ti(z) is the separation between two steps i and 
j, and i>2 = V20 (T — T t ) is the effective two step contact (short range) interaction. H int is 
the additional interaction, and two possible choices are 
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Hint = V 3 J J2 S ( r ij( Z )) S ( r ik(z))dz, (4b) 



H int = hJ2 j I r iiO) I dz. (4c) 

%<j 

The first form represents a three step contact repulsion while the last form represents a two 
step long range repulsive interaction. For T > T t , v 2 > and the noncrossing condition is 
ensured by taking the limit v 2 — > oo. This however is not required because it is known from 
RG that the repulsive case is described for d < 2 by a stable fixed point (FP) (see below). 

In a mean field treatment J2 ${ r ij) i n Eq. |4a] can be replaced by p 2 . Similarly the three 
body interaction would generate a p 3 term ||||, as also the r~ 2 repulsive interaction H5|Jl3| 
so that the mean field free energy is 

Fint{p, T) = f p + v 20 (T - T t )p 2 + cp 3 , (5) 

where c depends on h or v 3 as the case may be. This again gives (3 = 1, when Eq. |3| is used, 
though the physics behind this is completely different from that proposed is Ref. 0,0,0- 

Let us first consider the short range case, Eq. |4b[ A simple dimensional analysis shows 
that t>3 is marginal in d = 1. We introduce the dimensionless parameters u 2 = v 2 L 2 ~ d and 
u 3 = v 3 L 1 ~ d where L is an arbitrary length scale in the transverse direction. A renormal- 
ization procedure would take into account the effects of interactions at scales < L along the 
steps, changing the effective interaction felt at length scale L ("coarse graining"). The details 
can be found in Refs. fl4]|],|7[]. The running coupling constant for the two step interaction 
is known exactly [T^J^,^] and is given by 

Lj£ = (2 - d)u 2 - u 2 2 /2n. (6) 

The flow of u 2 is controlled, for d < 2, by the two fixed points u 2 = (unstable) and 
u 2 = 2n(2 — d) (stable). The unstable FP corresponds to the transition point for two chains 
and the tricritical point in the many chain case, while the stable FP describes the repulsive 
steps acting like fermions in d = 1. For u\ = 0, the RG equation for u 3 is given by H§] 

t du3 2 / 7 X 

l ~qI = -°3U 3 , (7) 

with only the fluctuation contribution in the higher order of u 3 at d = 1. [|TJ| Around the 
unstable fixed point for u 2 , for small deviations, the effective coupling is given by u 2 (L) ~ 
u 2 L for d = 1, so that the renormalized but not rescaled coupling constant is just v 2 . 
In contrast, the renormalized three step interaction at L ~ p" 1 gives v 3 (L) ~ v 3 /lnp. 
Substitution of these changes the p 3 term of the free energy of Eq. [| to p 3 / In p. The shape 
of the coexistence curve is then 

v 2 ~ p/ In p, i.e. p ~| T — T t \ In | T — T t \ . (8) 

We see that fluctuations produce a zero slope coexistence curve, though the coexistence 
exponent (3 is still 1, the mean field value [pi]] . For the high temperature phase, the system 
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is described by the stable FP u\ = 2 — d, and i^(L) ~ u\p yielding the famous p 3 term that 
produces the Kasteleyn-Pokrovsky-Talapov transition. [12|,19 



We can also predict the behaviour right at the tricritical point. With t> 2 — 0, the 
analogue of the Pokrovsky-Talapov transition would involve only the three body repulsive 
interaction. The relevant behaviour comes from the minimization of the free energy F = 
(fo~ p)p+V3p 3 / logp. Therefore, the step density at tricriticality behaves like p ~| p — fo \ 1 ^ 2 
(log | p — fo I) 1 / 2 with /o < 0. The exponent is the same Pokrovsky-Talapov one but with 
an additional log correction (which may be hard to detect). 

It is possible to have higher order multicritical points with just t> 2 and v m > involving 
an m step repulsive interaction (m > 3). An exponent of j3 = 1/2 can be recovered || for 
m = 4 in an mean field way because d — 1 is above the upper critical dimension of v m . Such 
a multicritical point requires v 3 = and with nonintersecting steps, it seems very unlikely 
that this will happen. 

Let us now come to the long range repulsion case, Eq. [|c[ LR interactions are special by 
virtue of their singular nature. A renormalization group transformation is analytic in nature 
and, therefore, can never generate a singular potential. A corollary of this is that the two 
body LR interaction does not get renormalized but affects the renormalization of the short 
range (non singular) pair potential. Such a renormalization is going to change the exponent 
of p in the v<i term of Eq. [5], and, therefore, the nature of the coexistence curve. The RG 
approach for this case is also available in the literature, and we quote the results |22|fl. 
Defining u{L) = a[v 2 (L)L 2 - d + h(L)A], g(L) = 2K d Ah(L), where A = 2 d - 1 n d / 2 T(d/2), K d 
is the surface area of a unit d— dimensional sphere, and a is a (system dependent) constant, 
the recursion relations from Ref || are 

L du/dL = —{u — u* s ){u — w*), and dg/dL = 0, (9) 

where u* s u = {2 - d ± [(d - 2) 2 + 4g} 1/2 }. Since g is marginal, and the FPs for u depend on 
the long range coupling g, we find a nonuniversal behaviour. The unstable FP u* u describes 
the tricritical point, so that linearizing around it, we can determine the effective coupling 
that goes in the free energy, provided g < 3/4. For two chains with g > 3/4, the binding 
transition is first order ||. We assume that the tricritical point, as an end point of the 
coexistence curve, has some critical nature. Therefore, g < 3/4. A straightforward analysis 
then gives Au{L ~ p" 1 ) ~| T - T t \ p 1 ^, where x = u* s - u* u = [(d - 2) 2 + 4s-] 1 / 2 . This 
gives, via Eq. |^ a coexistence curve p ~| T — T t p 3 , where, to repeat, x is a nonuniversal 
number. In order to achieve consistency with experiment, one requires x ~ 2, which in turn 
requires g ~ 3/4. The RG analysis of Ref |J, as already mentioned, also shows g = 3/4 
is a very special point, corresponding to an "upper critical dimension" case. Furthermore, 
for g = 3/4, with hard core repulsion, log corrections are expected, which are not captured 
in the simple RG analysis PJ23| . We conclude that if h of Eq. |4c] happens to be close to 



3/47T, then the coexistence exponent j3 will be close to 1/2, and if h = 3/47T, (3 = 1/2 with 
additional log corrections. 

In both the SR and LR cases, since the free energy is known, the surface stiffness can also 
be calculated. Following Ref. 0, we find that at a given T > T t , for the short range case, 
the surface stiffness approaches the free fermion value in a singular fashion ~ 1/ | ln6/d | 
with the misorientation while for the long range case, the free fermion value is reached from 
above in a 9 independent way. 
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To summarize, we considered two different scenarios both of which give a zero curvature 
coexistence curve. The purely short range interaction however only gives a log correction 
to the mean field exponent which seems to be far off from the experiment. The long range 
case predicts a nonuniversal value and the observed exponent seems to suggest that Si(113) 
surface with a miscut towards [001] is at or close to a very special point for the long range 
interaction. It is rather striking that the very first system that showed the tricritical point 
also corresponds to the very special point for the long range interaction. We are not sure 
whether it is just an accident or a general rule. If an accident, then other orientational phase 
diagrams should be studied experimentally to verify the claim of nonuniversality (and may 
be a simple verification of RG in statistical mechanics). If not an accident, we wonder why 
nature chooses to be at the threshold. 

I thank Sutapa Mukherji for many helpful discussions. The relevance of the long range 
case was suggested by M. Lassig whom I also thank for discussions at the initial stage of 
this work. This work is partially supported by DST SP/S2/M-17/92. 
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FIGURES 

FIG. 1. (a) Schematic diagram of the steps of equal height d. The misorientation is related 
to the density p, as tan^ = d/l = pd. (b) The steps viewed from above, (c) Schematic phase 
diagram. 
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